For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Hibbos1972 


THE UNIVERSITY OF ALBERTA 


A STUDY OF APERIODIC, SPATIALLY FINE VARYING TWO- 


DIMENSIONAL SOURCES IN ELECTROMAGNETIC INDUCTION 


by 
co) ROY DEAN HIBBS, JR. 
A THESIS 


SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN 
PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF 


MASTER OF SCIENCE 


Department of Physics 


Edmonton, Alberta 


Pails i. 


7 : —— 
Sn 
Bee Relay sel): eT Saas 0H lag" 


a” a5275% cat OE TRU ies Jecotcasere J 
F ah 
4 6 te - — 
ait Sa - : 


> 'See726904 


za] eat 2000 48 


UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE 


The undersvgneo certity that, they 
Faculty of Graduate Studies 
a thesis: entitled 2 STUDY GF 
SPATIALLY TIME VARYING 


ELECTROMAGNETIC. INDUCTION, 


thepetton £liahinent of the requirements co 
On, Master Of SCLEence:. 


SlUDIES AND RESEA 


TWO-DIMENSIONAL SOURCE 


RCH 


have read, and 
and Research 
APERIODIC, 


5S IN 


submitted by Roy Dean Hibbs, Jr. 


r the degree 


ABSTRACT 


Electromagnetic induction in the earth by a 
generalized two-dimensional source is considered. The 
global and local preblems are discussed in detail. A 
description of the two-dimensional E and H polarization 
cases is presented along with a review of the determinacy 


and the various methods used to solve the problen. 


A comprehensive discussion of the source is 
given and solutions for selected non-uniform sources are 


described. 


The electromagnetic induction by a symmetric 
ionospheric current distribution which is sinusoidally 
time dependent over a horizontally stratified earth is 
investigated and a method is presented for the calculation 
of perturbations of such fields by embedded 
inhomogeneities. This method is then extended to non- 
symmetric current distributions over layered and laterally 


inhomogeneous subsurfaces. 


An aperiodic source is obtained by summing the 
Single frequency solutions for layered subsurfaces ina 
Fourier transformation. By superposition of the 
electromagnetic fields of separate aperiodic, non-uniform 
current sources which vary independently with time a 


source which exhibits spatial variation with time is 
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1. INTRODUCTION 


1.1 General 


The inductive response of the electrically 
conducting earth's interior to external gecmagnetic 
Variations can be inferred from the observed surface 
fieids. Currents induced in the earth by external 
oscillating electromagnetic fields considerably modify the 
surface field values. These external sources consist of 
current systems of finite extent flowing in the ionosphere 
or above. In general the geomagnetic variations produced 


by these current systems are transient in time. 


In this thesis the global and local problems in 
electromagnetic induction are discussed and the general 
two-dimensional local perturbation problem is considered 
together with various methods for its solution. The 
source effects due to non-uniform source current 
distributions are then investigated and solutions for 
selected non-uniform sources over a layered earth are 
presented. A method for solution over a laterally non- 
uniform earth is described. Symmetric and non-symmetric 
ionospheric current distributions are considered and 
solutions for such current distributions which vary 
Sinusoidally with time over a uniformly layered earth as 


well as over a laterally inhomogeneous earth are computed 
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and compared with uniform source results. The 
mathematical formulation of the boundary conditions 
necessary to obtain these solutions is given. Aperiodic 
sources are investigated by the summation of the solutions 
for different frequencies through the use of the Fourier 
transform. By the superposition of solutions for 
different current distributions which do not vary 
Sinusoidally in time, a source which changes spatially 
with time can be obtained. The magnetic field components 
of such a source over a layered earth of finite 
conductivity are then computed and profiles of the 


components obtained. 


1.2 Global and Local Problems 


The mathematical problems encountered in the 
study of earth currents are of two types. The first type 
involves inducing fields of global dimensions with 
averaged conductivities for the earth as a whole; these 
are known as global problems. The second type may or may 
not involve inducing fields of global dimensions but the 
conductivity values are definitely local values; these are 
known as local problems. Following the discussion by 


Price (1964) these problems are further considered. 
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1.2.a Global Problems 


In the first group of problems, since we are 
interested in large regions having dimensions comparable 
to those of the earth, it is necessary to use spherical 
polar coordinates. The electrical conductivity v is then 
treated as a smoothed function of the coordinates (fr,0,0) 
at any point within the earth. This function will not 
take into account any local variations in conductivity 
(which may be considerable) but only large scale 
variations of suitably defined average v for any region. 
For this problem the currents induced in any earth model 
will depend on the nature and distribution of the inducing 
field. The first step in solving such a problem is to 
express the inducing field in terms of spherical harmonics 
(assumed given). The induction effect of each harmonic is 
then studied separately and the result obtained from the 
summation of all the harmonics. The basic theory and 
method of solution for this type of problem is straight- 
forward, though the calculation may be difficult and 


laborious. 


One global problem which has been extensively 
studied is that of the induction, by a varying magnetic 
field, of electric currents in a spherical conductor 


in which the conductivity is a function of r, the distance 
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from the center of the earth. Lamb (1883) treated the 
case where v is uniform, and Lamb's solution was used by 
Schuster (1889) in his discussion of the internal and 
external parts of the magnetic daily variations. Chapman 
(1919) and Chapman and Whitehead (1922) made further 
applications of Lamb's solution. Price (1930, 1931) 
extended Lamb's solution to aperiodic fields and Chapman 
and Price (1930) then used this extension to discuss the 
induced part of the field of magnetic storms. The more 
general case when v is any function of r was studied by 
Lahiri and Price (1939). They developed both formulae and 
methods of calculation for both periodic and 

aperiodic fields when dealing with any distribution of 
conductivity within the sphere r=a that can be expressed 


in the form: 


v=v9 (qa/r)” for r<qa<a and v=0 for r>ga (1.1) 


where v9, gq and m are constants. 


Lahiri and Price (1939) applied their results to 
investigate the earth's conductivity to considerable 
depths and many applications of their formulae and methods 
have since been made by other investigators (e.g. 


Rikitake, 1950). 
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Apparently there have been no successful 
attempts to obtain solutions (analytical or numerical) for 
the induction of currents in a spherical conductor, when 
the conductivity is not only a function of r but also 
depends upon --or @ or both. It seems reasonable that, 
except within the outermost 700 km. of the earth, the 
variation of the conductivity with r is far more important 


then any dependence upon + or @g. 


The variation of v with @e© and @ within the crust 
is of importance in both global and local problems. For 
example, the conductivity of the oceans is several orders 
of magnitude greater than that of the continents. Non- 
uniform spherical shells are useful in modeling the giobal 
effects of currents induced in the crust. Price (1949) 
has given the basic theory of electromagnetic induction in 
non-uniform thin sheets and shells and applications have 


been made by Ashour (1950), Rikitake (1960), and others. 


For very slowly varying magnetic fields it is 
not possible to treat the crust as a thin shell by itself, 
because of the mutual induction between the crust and 
conducting mantle below. This point was illustrated by 
Rikitake (1961) in his discussion of the currents induced 


in the oceans by daily variations. 
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1.2.b Local Problems 


In the local problems the sphericity of the 
earth is ignored, and the earth is treated as a semi- 
infinite plane body having a non-uniform distribution of 
conductivity, where our attention is confined to some 
limited portion of the conductor. However, the inducing 
field is usually of large dimensions and effectively 
uniform over the region of interest. This has sometimes 
been taken to imply that the actual distribution of the 
inducing field can te ignored when calculating the induced 
currents and fields in the region considered. fThis 
assumption is implicit in the theory of the 
magnetotelluric method of geophysical prospecting given by 
Cagniard (1953). The limitations of results based on this 


theory are discussed in section 2. 


Price (1964), in a comprehensive note on the 
interpretation of magnetic variations and magnetotelluric 
data, points out that a knowledge of the inducing field 
and surrounding ground conductivity in the immediate 
neighborhood of some particular station is not usually 
sufficient to determine the strength of the induced 


currents in that neighborhood. The induced currents are 
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affected by the distribution of the entire inducing field 
and by the average properties of the conductor over a 
region of corresponding dimensicns. It is true that the 
entire induced current system will also contribute to the 
field at a particular station, but the largest 
contribution comes from induced currents in the immediate 
neighborhood. The important point is that not only the 
local properties, but the properties of the conductor as a 
whole determine the strength of the induced currents in 


the region of interest. 


Price (1964) also points out that to approach 
the problem of the interpretation of local patterns of 
geomagnetic variations in an objective way, the variation 
fields should be separated into parts of external and 
internal origin and these parts should be examined 
separately. The anomalous patterns in the fields of the 
induced currents are produced by non-uniform distributions 
of conductivity in the earth, whereas the external 
inducing field is of a relatively simple form and varies 
only slighty over the anomalous area except near the 
equatorial and auroral electrojets. There seems to be no 
reason to suggest that the ionosphere has anomalous locai 
properties, even though a small effect on the (non- 


isotropic) conductivity of the ionosphere might be 
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produced by a strong perturbation of the earth's magnetic 
field. Price (1964) suggests that even if a considerable 
disturbance is caused in the ionosphere by such local 
anomalies it is doubtful that the surface values of the 
magnetic field would be greatly affected. Also, if the 
external field is found to reflect to a small extent the 
anomalous magnetic field distribution, this can be 
attributed to the general mutual induction between the 


earth and the ionosfhere. 


Once the external field and its induced current 
system is known, the local problem reduces to examining 
the local redistribution of a given average system of 
induced currents caused by local inhomogeneities of 
conductivity. However, it is not necessarily a simple 
problem of redistribution of steady current flow, since 
the currents are time varying and therefore a skin effect 


will be encountered. 


The mathematical problems which must be 
considered in this regard are those relating to the 
disturbance of skin effect distributions of current rather 
than actual electromagnetic induction problems. This type 
of mathematical problem is described by Price (1964) as 


follows: 
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"Using cartesian coordinates (x,y,Z) with the z 
axis vertically downwards, a non-uniform conductor 
occupies the half-space z>0. Near the origin the 
conductivity is a function (not necessarily continuous) of 
(X,Y,Z), but at great distances from the origin it is a 
function of z only. A given alternating e.m.f. impels 
currents near the surface of the conductor. [The problem 
is] to determine the distribution and surface field of 


these currents." 


The electric field arises from a varying 
external magnetic field. This electric field is 
controlled not only by the external source but also by the 
properties and geometry of the conductor at great 


distances. 


Problems of this sort can be solved for two- 
dimensional cases but when an attempt is made to reduce 
any problem involving electromagnetic induction of 
currents to a two-dimensional problem by supposing that 
conditions obtained in a limited region extend to infinity 
in two opposite directions, it is necessary to closely 
examine the conditions at infinity and ensure that they 


are acceptable. 


The problems of electromagnetic induction dealt 
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with in this thesis are confined to two-dimensional 
problems with regard to both conductive geometry and field 


adistribution. 


1.3 The Two-Dimensional Probiem 


1.3.a E and H Polarizations 


The coordinate system is shown in Fig. 1.1. The 
region -|h|<z<0 is a free space region with conductivity 
equal to zero. A conductor occupies the half-space z>0 
and it is assumed that a current source of some form is 


located at or above z=-{h|. 


For an oScillating field with time dependence 
exp (iwt) where w is the angular frequency and the period 
2W/w is sufficiently long so that displacement currents 


may be neglected, Maxwell's equations become: 


CURT RH y.e (1.1) 


CURL E = -iwuH (1.2) 


where v is the conductivity appropriate for each region, u 
is the magnetic permeability of free space and the time 


factor exp(iwt) is understood in all field quantities. 
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In the two-dimensional problem the fields do not 
vary in one direction, in this case the x-direction. That 
is all guantities are independent of x, and Maxwell's 


equations reduce to 


dgHZ/Jy - dHY/Q2 = VEX, (4.3) 
Pixs 2o=.vey; (1.4) 
-)Hx/Jy = vEz, (1.5) 
aEZ/aY.- ag HY/I2 = iwuHkx, (1.6) 
g EX/)Z = -iwuHy, (1.7) 
—JEX/)Y =. -iwuHlz. (1.8) 


These equations are such that only Ex, Hy and Hz 
aresSinvolved in (1.3, 1.7 and 1.8) and only Hx, Ey and Ez 
are involved in (1.4, 1.5 and 1.6). These two separate 
sets of equations can therefore be solved independently. 


Combining (1.3), (1.7) and (1.8) we obtain 


a°EX/Jy* + J2ExX/) z2 = iwuvEx (129) 


which is the equation to be solved in all regions for the 
E-polarization case (the case in which the electric field 


is parallel to any conductivity discontinuities) . 


Combining (1.4), (1.5) and (1.6) a Similar 


equation is obtained for the H-polarization case (the case 
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in which the magnetic field is parallel to any 


conductivity discontinuities) : 


ag *HxfAvy? + g2Hxfyze = iwuvHx. (1.10) 


With the proper values of v inserted for each 
conductive region, and with the proper boundary 
conditions, these equations can be solved for the separate 
cases to obtain the field distributions. For the E- 
polarization, Ex is determined from the solution of 
(1.9), and then (1.7) and (1.8) may be used to calculate 
the appropriate components. Also, for the H-polarization, 
Hx may be determined from (1.10), and equations (1.4) and 
(1.5) may be used to obtain the appropriate electric field 


components. 


1.3.b Determinacy of the Two-Dimensional Problem 


As mentioned in section 1.2.b not all two- 
dimensional problems are determinate. It is important, 
then, to examine the conditions for determinacy for a 
periodic inducing field over a semi-infinite uniform 
conductor with a plane boundary occupying the half-space 


z>0. Following Price (1950), we let Ex be of the forn 
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EX = Z(z).Y(y).exp(iwt), (1.11) 


and substitution of this expression for Ex into equation 


(1.9) gives: 


O27 T= (I7 Sy). (iwuvZ - 9 227) z=) .Y Clint) 


fOorez eo ,-and 


d?7Y/Jy? = -(1/2).- (9 22/9 z?).Y (leal3) 


for z<0, where the arguments for Y and Z have now been 


dropped. 


Since the coefficient of Y on the right of 
(1.12) and (1.13) is independent of x and y it must be a 
constant, -s?, where s is real, positive and non-zero 
such that 
a7Y/ y* = -s#Y. (1.14) 


Therefore, in the region z>0 


d2Z/)z2 - (iwuvt+ts2)Z2 = 0 (ia) 
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and 


Z(z) = A.exp(-pz) + B.exp (pz) 


where p2=iwuv+s2 and A and B are constants. 


In the region z<0, 


)2Z/)z2 - s2Z = 0 


and 


Z(z) = a.exp(-sz) + b.exp(sz), 


where a and b are constants. 


The corresponding E field is then 


ExX = {A.exp(-pz) + B.exp(pz)}.Y 


for z>0,. and 


Ex = {a.exp(-sz) + b.exp(sz)}.Y 


ZOD 2404 


As z becomes large, Ex must tend to 
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B must be zero. At the boundary z=0, Ex must he 


continuous, and so 


Ay= artob: (¥.24) 


The corresponding magnetic fields are obtained from (1.7) 


anagmGUac Lor 20 : 


Hy = (i/wu).(gExf~z) = -(ipA/wu).exp(-pz).Y (1.22) 
Hz = -(i/wu). QEx/)y) = - (iA/wu).exp (-pz) 
-QYAY); (1.23) 


and Lor z<0; 


Hy -~(is/wu). f{a.exp(-sz) - b.exp(sz)}.Y (1.24) 


HZ -(i/wu). fa.exp(-sz) + b.exp(sz)}.@Y/y). (1.25) 


The tangential component of the magnetic field 
must be continuous across the boundary z=0, which leads to 


the relation 


pA = s(a-b). (1.26) 


The scalar potential & of the total magnetic 


field in the region -|h|<z<0 satisfies Laplace's equation 
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& = —{C.exp (Ssz) ) + Dsexp(sz).}.P (y) 


where P(y) satisfies the equation g2P(y) J y2 
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Since the conductivity outside the conductor is 


Zero 


and Hy and Hz can be found from the relation 


H = -grad & 
as: 
Hy = {C.exp(-sz) + D.exp(sz)}.-(dP(y) AY) 


HZ =S. {C. exp (—SZ) = Dw exp (SZ) } oP (y).« 
Equating (1.24) and (1.30) gives: 


-isa/wu = C isb/wu = D Y(y) = gdP(y)AYy- 


COMDANLNG) (e21).,.90t.20)) and (1.32) (gives; 


(7528) 


(130) 


(Ves) 


(1.32) 
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D= {(P-S)/(PtS)}-C, (1.33) 
de= et (Dts) / (prs) }. (iwu7s) =D, (1.34) 
b= -{(P-S)/(pts)}.-a, (1-35) 
and 

Ree —2S)D/ (pS). (1.36) 


Therefore by specifying C all the other 
quantities will be determined. This anaylsis, while 
confined to uniform subsurfaces, can be extended to multi- 


layered cases. 


The part of & in (1.27) which involves exp (-sz) 
corresponds to the inducing field of the current sources 
and the part which involves exp(sz) corresponds to the 
field of the induced currents. It can be assumed that the 
inducing field potential is specified by the current 
distribution so that C.exp(-sz).P(y) is known. From the 
preceding discussion it can be seen that all other fields 
and potentials are determined from the inducing field. 
Therefore, the problem is determinate for s positive and 
non-zero. For s=0 the problem becomes indeterminate in 
general as explained by Price (1950). The H-polarization 


Can be shown to be determinate in the same manner. 
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1.4 Methods for Solving the Two-Dimensional Problem 


In his classic paper, Price (1950) considered 
electromagnetic induction in a semi-infinite conductor 
with a plane boundary for any inducing field. Cagniard 
(1953), in his derivation of the well known expression for 
apparent resistivity, considered a horizontally stratified 
earth and a uniform source. Wait (1954) and Price (1962) 
discussed Cagniard's results in terms cof a source of 
finite dimensions. However, until the last decade, little 
attention has been paid to the problem of a vertical fault 
in the conductive region which represents the earth. It 
is well known that there are limited regions in various 
parts of the earth where temporal changes in the 
geomagnetic field at relatively near stations show 
considerable differences in form and amplitude. These 
differences remain consistent in character for the 
fluctuations of a given frequency, and forma 
definite pattern for the area studied. fThey are 
associated with lateral variations in conductivity below 
the earth's surface. It is important to understand the 
manner in which lateral variations affect the surface 


field in order to interpret such variations. 


D'Erceville and Kunetz (1962) obtained a 


solution for two media of different conductivities in 
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contact along a vertical plane overlaying a horizontal 
basement that is either infinitely conducting, infinitely 
resistive or at infinite depth when the magnetic field is 
everywhere parallel to the strike of the fault, (H- 
polarization). Rankin (1962) extended this work in his 


study of "conductive" and "resistive" dikes. 


Weaver (1963) considered a half-space conductor 
which consisted of two quarter-spaces of different 
conductivity with a plane vertical fault of infinite 
depth. He considered both the H-polarization case and the 
E-polarization case using an analytical technigue. His H- 
polarization case agreed well with that of D'Erceville and 
Kunetz. However, in order for him to obtain a solution 
for the E-polarization case he found it necessary to use 
the approximate boundary condition that the tangential 
magnetic field is constant along the surface of the 
conductor. Recently, Mann (1970) has proposed a 
perturbation technique to consider the original Weaver 
(1963) problem, and has shown that Weaver's original 
solution is just a first approximation to the field inside 
the conductor for normally incident waves. Weaver and 
Thomson (1972) have applied this technique and have been 
able to avoid the approximate boundary ccndition of the 


constancy of the horizontal magnetic field along the 
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surface and have obtained approximate solutions. 


In the meantime several numerical methods have 
been used to obtain solutions for problems with vertical 
discontinuities for both Pep lacie aed and E-polarization 
cases. Dulaney and Madden (1962) were the first to apply 
the transmission line analogy over a two-dimensional mesh 
to such problems as these. This method has been used by 
many authors (including Madden and Thompson, 1965; 

Madden and -Switt, 1969; Swift; 196/,1971<2 Wright, 
1969,1970; and others). The finite element method, which 
uses the principle that electromagnetic fields behave in 
Such a way aS to minimize the energy, has recently been 
applied to induction problems by Coggon (1971), Ryu (1972) 
and Reddy and Rankin (1972). Finite difference techniques 
have been applied by Neves (1957), Latka (1966), Patrick 
and Bostick (1969) and others. Jones and Price 
(1969,1970,1971a,b) have employed the finite difference 
technique for studies of various two-dimensional 
conductivity distributions, and Jones and Pascoe (1971) 
and Pascoe and Jones (1972) have given a general computer 
program for the solution of the local perturbation 


problen. 


Jones and Price (1969,1970) considered the half- 
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Space conductor as the limit of a spherical conductor as 
the radius becomes infinite. They then carefully 
developed the boundary conditions and were able to obtain 
solutions for both the H-polarization and E-polarization 
cases. Jones and Price (1970) studied the field 
distribution within the whole two-dimensional region of 
interest, as well as the surface values of the various 
components, including the Cagniard (1953) apparent 
resistivity. In more recent work, Jones and Price 
(1971a,b) and Jones (1971a,b) as well as Hibbs and Jones 
(1972a) have calculated apparent resistivity vaiues along 


the surface of the conductor. 


It should be noted that in all the numerical 
methods mentioned above the assumption has been made that 
the inducing field is uniform. A non-uniform source in 
the form of a line current has been considered by Wait 
(1962), Dosso and Jacobs (1968) and others. Schmucker 
(1971) has considered non-uniform sources over a laterally 
inhomogeneous earth in which the lateral changes in 
conductivity are confined to a limited depth range. 
Hermance and Peltier (1970) studied the magnetotelluric 
field of a line current over a layered subsurface, and 
later (Peltier and Hermance, (1971)) extended their 


approach to include symmetric current distributions, such 
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as a Gaussian current intensity source. 


One aspect of the work of this thesis has been 
to combine the method of Peltier and Hermance (1971) with 
the iterative finite difference rernad of Jones and Price 
(1970) to investigate embedded conductivity 
inhomogeneities for symmetric non-uniform sources, as well 
as for non-symmetric current source distributions. Two 
papers (Hibbs and Jones, 1972b,c) describe this work and 


it is further discussed in the following sections. 


It should be emphasized here that the E- 
polarization case is readily adaptable to such an 
approach, whereas the H-polarization case is not. 
Consequently, work of this thesis is primarily concerned 


with the two-dimensional E-polarization case. 
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2. , SOURCE DISCUSSION 


Cagniard (1953), in his classic paper presenting 
the basic ideas underlying a prospecting method using 
magnetotelluric techniques, assumed that natural 
oscillating electromagnetic fields near the earth's 
surface are uniform. Consequently, the amplitude and 
phase relationships of the horizontal components of E and 
H over a layered conductor depended upon the electrial 
conductivities of the strata and the period of the 
oscillating field. The apparent resistivity in ohm-m is 


given by 
p = 0.2T(E/H) 2 (2.1) 


where T is the period in seconds, and E (in 
millivolts/kn.) and H (in gammas) are mutually 
perpendicular electric and magnetic field components at 
the surface of the earth. Cagniard (1953) presented 
Master curves of apparent resistivity for various layered 
conductors and a method of interpretation of experimental 


magnetotelluric data from the master curves was given. 


In his original paper Cagniard (1953) did not 
explain the origin of the earth currents but only required 
that they be uniform. Later, in correspondence with Wait 


(which is reproduced following the paper by Wait, 1954) 
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Cagniard stated that the magnetotelluric perturbations are 
caused by vast systems of ionospheric currents whose 


dimensions are on a global scale. 


Since the depth of penetration of induced 
currents depends upon the period, by using one observing 
station the magnetotelluric method is able to indicate the 
conductivity distribution with depth (over a suitable 
range of frequencies) and a conductivity profile can be 
deduced. However, if the source field is non-uniform, the 
influence of such non-uniform source field distributions 
should be studied in relation to the magnetotelluric 
method. It is apparent that over some regions of the 
earth, for example near the magnetic equator (Forbush and 
Casaverde, 1961) and near the auroral oval (Kisabeth and 


Rostoker, 1971), the inducing field is not always uniform. 


Aside from the discussion by Wait (1954), little 
attention has been paid to Cagniard's fundamental 
assumption that the electromagnetic field is uniform over 
any horizontal plane. Wait (1954) showed that if this 
condition is not satisfied Cagniard's calculation of the 


quantity E/H is not exact. 


From calculated corrections to this quantity 


involving the second-order spatial derivative, Wait 
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concluded that the corrections would be necessary if the 
magnetic field changed appreciably in a distance of 35 
km. for periods greater then 10 seconds and ground 
conductivity of the order of 10-3 ohm/m. Wait also 
observed that if the source field was generated by 
ionospheric currents flowing at heights of about 100 

km. the corrections would be important for many of the 
geomagnetic oscillations used in magnetotelluric methods. 
This is particularly true of geomagnetic fluctuations 
observed in connection with the fields of auroral and 


equatorial electrojets. 


As mentioned before, Cagniard argued in his 
reply to Wait that the magnetotelluric perturbations are 
produced by vast current systems whose electric fields 
have global dimensions, and therefore his equations could 
be applied in most cases. This conclusion is fairly 
reasonable for the source fields considered by Wait. 
However, Wait himself made a simplifying assumption in 
that he considered a semi-infinite conductor of uniform 
conductivity. For a conductor whose conductivity varies 
with depth it is found (Price, 1962) that the dimensions 
of the source field cannot be ignored even if the field is 
of global dimensions and the skin depth of the inducing 


field is moderate. 
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Price (1962) outlined the general theory of 
magnetoteliuric methods including source field 
considerations and applied this theory to a simple case to 
show the importance of the spatial dimensions of the 
source. He estimated the spatial wave number s of the 
source field and concluded that the least value of s can 
be obtained by equating 2W/s to the circumference of the 
earth, which corresponds to the wavelength of a field 
representing a spherical harmonic of the first order, 
which gives s=1.57x10-9 cm+1. Most geomagnetic fields 
could contain spherical harmonics of a higher order than 
this, and a more representative value would be s=1.0x10-8 
cm-1, More localized fields such as those of an 
ionospheric electrojet contain spherical harmonics of an 
even higher crder. Price suggests that a representative 
value of s could be found by equating 2W/s to four times 
the height (approximately 100 km.) of the ionospheric 
currents. The maximum value of s is thus approximately 
1.57x10-7 cm-1. Thus the values of s of most interest in 
magnetotelluric investigations lie between 1.57x10-9 cm-! 


and W.5/7x10—-7 cm-+. 


Price (1962) then showed, that for certain 


distributions of subsurface conductivity, the quantity 
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Ex/Hy can be considerably affected by the value of s. One 
Simple case which he chose was that in which the conductor 
has a constant conductivity v down to a certain depth D 

and zero conductivity below. The gquantity Ex/Hy was found 


to be 


ExX/Hy = (iw/p). {{ pts+(p-s).exp(-2pD) J/[ pts- (p-s) 


exp (-2pD) J} (2.2) 


where p*®=s@+tiwuv. 

From this it is seen that the ratio Ex/Hy depends upon s, 
both explicitly on the right hand side of (2.2) and in the 
quantity p as well. When s=0, which corresfonds to a 


uniform source field, 

Ex/Hy = (iw/pe). {{ 1texp(-2peD) J/[ 1-exp (-2peD) J} (2.3) 
where pe=/iwuv. 

When D is infinite, that is for a semi-infinite conductor, 
equation (2.2) reduces to 


EX/HY = iW/p (2.4) 


and so for this case there is a dependence on s. 
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With a first layer resistivity of 1000 ohm-m at 
a period of 100 seconds, the amplitude of |{Ex/iwHy| was 
found by Price (1962) to vary greatly with spatial wave 
number s which itself varies over the range 1.57x10-9 cm-1 
to 1.57x10-7 cm-1. This was particulary evident when the 
thickness of the first layer of the conductor was less 
then 10 km. A complete description of this phenomenon is 


given by Price (1962). 


Although the foregoing calculations were made 
for a simple distribution of conductivity, they show that 
the dimensions of the source field can be an important 
influence on the quantity |Ex/iwHy|]. This is true not 
only for local fields but also for fields with global 
dimensions. It’ should be pointed out also that the above 
discussion deals with sources of a single wave number 
whereas sources in general contain many wave numbers. The 
sources investigated in the remainder of this thesis will 


deal with the latter more general type. 
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3. SOLUTIONS FOR NON-UNIFORM SOURCES 


3.1 Two-Dimensional Solutions for Non-Uniform Sources over 


Layered Media 


3. tea Laine Current Source 


Hermance and Peltier (1970) investigated the 
magnetotelluric fields of a line current over a layered, 
two-dimensional earth. The line current flowed ata 
height of z=-{h{ parallel to the x axis of Fig 3.1 with a 
Sinusoidal time dependence exp(iwt). The E-polarization 
case is the applicable two-dimensional problem in this 
regard, and for this case the equation to be solved is 


given by (1.9). By the method of separation of variables, 


solutions of the form Ex(y,z) = Y(y)Z(z) are obtained, 
where 

Y(y) = C(s).cos(sy) + D(s).sSin(sy), (3-1) 
Z(z) = An(P).exp(-—Pz) + Bn(P) -exp (Pz), (82) 


and P*=s@+iwu.Vn. 
The values of s and Vn are the spatial wave number and the 


homogeneous conductivity of the nth layer respectively. 


Due to the position of the line current the 
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electric £veldvis®anveven function of y. ‘Therefore only 


the cosine term of (3.1) need be considered. 


In the space above the conductor the 
conductivity is zero and PLTSs If 4, (P)=1, then C(s) is a 
frequency dependent parameter which characterizes the 
source. ep ed can be thought of as a reflection 
coefficient which represents the contributions to the 
total field from subsurface layers (z>0). From the above 


discussion the elementary electric field in the first 


layer is 


Ex(y,zZ) = {exp(-sz) uy) (P) .exp(Sz)}.C(s) 


«COS (SY) « (3. 3) 


The actual field components are synthesized from 
the elementary fields by the use of the Fourier integral: 
ve 
EX (y, Z) = f texp (-sz) + B, (P) «exp (sz) }-C (s) 
Ce] 
-cos(sy) ds. (3.4) 
By taking the curl of (3.4) and applying 


Maxwell's equations the magnetic field components are 


found: 
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-€(s).cos(sy) ds, (325) 
Hz(y,Z) = (i/wu) . "s. {exp (-s2) + Bye Bi ene (Sot) 
0 


~C(s)issin (sy) ds. (3.6) 


If the contribution frem the subsurface layers is 
neglected by setting a) = 0, then a line current at 


z°=-h has a horizontal magnetic field given by: 
Hy¥@ee~ (2°72). (2-29) /(y*+ (2-29) 4) (3.7) 


where Ie is the line current intensity. 


Setting (3.5) equal to (3.7) gives 


reg 


~(i/wu) . { s.exp (-sz) .C(s) «cos (sy) ds 


er le/2n ye. (za2o)7 (yet (2-29) 4). (3.8) 
from Laplace transform tables (Erdelyi et al., 1961) 


exe (-L2-20 Js) «cos (sy) ds = (z-z9)/(y2+[z-z° ]#), (3.9) 


re) 


and combining (3.8) and (3.9) and solving for C(s) gives: 


C(s) = -(iwule/2T). (exp(z%s) /s). (3.10) 


The remaining unknown Bn(P) can be determined by 
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applying the boundary conditions that the horizontal 
electric and magnetic fields are continuous across the 
interfaces between the conductive layers. For a two 


layered earth B (P) must be zero in order that Ey3 does 


not become infinite with depth. These conditions give: 
eye = (R-1)/(R+1), (3.11) 
where 

es S(Q+1)/{P,. (Q-1)}, (3.12) 
and 

Rami Ane P,)}.exp(2P, d), (3.13) 


where d is the thickness of the first layer of the 


conductor. 


The relations for the electromagnetic field 
components for the first layer (3.4), (3.5), and (3.6) can 
then be numerically integrated on a digital computer using 
Filon's Method (Tranter 1956). The observed value of any 
field component is then obtained by taking the real part 


of the product of the field component and the time 


: - - 
. téi 4 sito ony 


- ' ene tense yi 
re ee | 


ms, : 


ae 


i ou 
| oi 

gaoh vq s8ef “Mh em 2h. C30" it sac (¢). 2 seas b 
toveg sucess dimen é4odtl | .aegeh EW etleitat es0os nae 


pues ah 


“pete 467 Jensen uve 


(rt .8) . yt rent —a) = 4 
an 
; 


7 
- =f 


(on st) | oettenp. ays beohae a 
*s. ca 7 


_ 


: i) 
. oj 4 = (a4 +.8)} = 
(€ fet) thy nPaRinwk | oie +g 
- 


a* %*>. geval 208 249 io egsatoda? ear 


= 


Blais yitenpemvisoels edt 103 stotzeden sif 

aes (0.6) ‘baw 5 (20%) 9 (O00) tayes Seige iz 

cAfay sovoqeos stip’ 5. go peep 896% ” 

tas. lo onlay wr deste ont a 520 rs 14d t5% 
dale lewt eds pains a be es ie 


33 


dependence exp(iwt). 


By varying the height of the line current, 
Hermance and Peltier (1970) concluded (with some 
reservations) that for an ‘average continental! crust 
Cagniard's magnetotelluric method is applicable if the 
electrojet can be represented by a line current ata 
height greater than 500 km., or by a broad current 
distribution at lower altitude. However, Hermance and 
Peltier (1970) felt that more extensive research into the 
validity of these restrictions would have to be 
undertaken, and concurred with Price (1962), Wait (1962), 
Ward and Morrison (1966), Rikitake (1966) and Madden and 
Swift (1969) that the current source structure could not 
in general ke ignored. Magnetotelluric measurements over 
tectonically active regions where highly conductive 
regions occur at depths of less than 25 km. were found to 
be much less susceptible to source effects, and therefore 
in these regions Cagniard's method could be applied 
(Hermance and Peltier, 1970). These observations are 
restricted to layered earth models and the source effect 
could be enhanced by lateral inhomogeneities in the 


conductivity configuration of the model. 


3.1.b Gaussian Current Source 
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AS an extension of their work on line current 
sources Feltier and Hermance (1971) considered a laterally 
diffuse current distribution. The electrojet was 
simulated by a continuum of weighted line current elements 
whose fields where superimposed to obtain the total field. 
An infinitesimally thin current sheet at a height of 
z=-{h{=-110 km. represented the source. The distribution 
of the current intensity was taken to be Gaussian and 
characterized by a standard deviation of 240 km. and the 


time dependence was sinusoidal. 


Since the current distribution is symmetric with 
respect to the origin, the electric field is again an even 
function of y. Following the same discussicn as 
previously, which led to (3.3), the elementary electric 


field in the first layer is 


Ex(y,Z) = {exp(-sz) + Bel eget (Sale t=) 


ecos (Sy). (3.14) 


B (P) is defined by (3.11) through (3.13) for a two 
layered earth and C(s) is obtained for a line current as 


before: 


C(s) =-{(iwule) /2™}. (exp (zs) /s). (S815) 
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The Gaussian current distribution can be 


represencted@asra = tunction@of Fy 
PAY) ls .expi(=y</ 2k=) (20716) 


where Ie and k are the current intensity directly above 
the origin and the standard deviation respectively. If it 
is assumed that the source is a composite of elemental 
line sources with an intensity distribution described by 
(3.16), then the horizontal magnetic field intensity due 


to an elemental line source at (y°,z°) can be written as: 


Hy?(element) = -(Ie/2W) .exp(-[y® j?/2k2). (z-z°) 


PRAZEE HEC aC ) eis (3.17) 


The total horizontal magnetic field is obtained by 


integrating over all values of y® and is 


co 
Hy!{total) = -(1e/2m). f {exp (“Ly ]2/2k2) . (2-29) 
fe) 


(ATA VIPERA Of Pe AN Be (3. 18) 


The Gaussian source constant G(s) is obtained in the same 


manner as C(s): 
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00 
-(i/wu) « fs.exp (~s2) «6 (s) «cos (sy) ds = -(Ie/2W). 
f) 


oo 
Jtexp (Ly? 12/2k2) . (2-29) 


Zip(ee2 wes (yay°) =p edy oe. (3.19) 
If in (3.9) y is replaced by (y-y®) and both sides are 


multipied by exp(-[y9]?/2k2) and the expressions are 


integrated over all y°%, then 


@ 00 
-(i/wuy ffs. exp (-sz) .C(s) «cos (sl ¥-¥°D 
( 
exp (-[y° ]2/2k2) ds dy® = -(Ie/2m) 


oo 


._f fexp (-Ly9]2/2k2) . (2-29) 


~o 


Feo EY Oe a ae yee (3.20) 
and comparison of (3.19) and (3.20) gives 


fe 2) 
c(s) {exp (-Ly° ]2/2k2) .cos (sl y-¥° ) dy® = G(s) 
— 


-cos(sy). (32211) 
By solving for G(s) we have 
G(s) =/2Wk.exp(-s2k2/2) .C(s), (3.22) 
and therefore 


G(s) = -(ilewuk// 2W.s) .exp(sz°-s2k2/2) . (3.23) 
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The electromagnetic field components for the 
first layer are obtained by substituting G(s) for C(s) in 


(3.4) 7. (2. Sjmande ere) % 


Again through numerical integration on a digital 
computer by Filon's method the actual field components may 
be obtained. The observed values of the field components 


are obtained as in 3.1.a. 


From the results of this work Peltier and 
Hermance (1971) concluded that the most pronounced source 
effect occurs over continental areas where layers of low 
resistivity overlay those of high resistivity. Over this 
type of conductor it was concluded that the conventional 
Cagniard method might lead to serious errors in 
interpretation. Also, it is possible that the presence of 
lateral inhomogeneities in the conductivity structure 
could enhance the effect of the non-uniform source. 
However, they found that the source effect for all two 
layer models in which a high resistivity surface layer 
overlies a lower resistivity layer is diminished for the 


Gaussian source distribution. 


Ae Two-Dimensional Solutions for Non-Uniform Sources 


over Laterally Inhomogeneous Media 
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3.2.a Iterative Solutions 


The primary limitation of the method of solution 
in the preceding sections 3.1.a and a: is that the 
conducting medium must be laterally homogeneous. While 
the solution for only two specific source configurations 
were obtained, other symmetric source configurations can 
be calculated and the method can be extended to non- 
Symmetric sources as will be shown in section 4.1. The 
solution for laterally inhomogeneous conductivity 
structures with finite conductivities can best be obtained 
through the use of the finite difference method. The non- 
uniformity of the source is reflected in the boundary 
values supplied to the finite difference mesh. A complete 
discussion of the boundary conditions will be given in 


Section S.2<D. 


Jones and Pascoe (1971) presented an iterative 
finite difference method for dealing with an embedded two- 
dimensional anomaly of arbitrary shape. This work was an 
extension to the method developed by Jones and Price 
(1970, 1971a,b) who considered several particular 
laterally inhomogeneous conducting regions. The method 


involves the solution of the appropriate finite difference 
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equations over a mesh of grid points by the Gauss-Seidel 
iterative technique (Smith, 1969). The general model is 
illustrated in Fig. 3.2 along with the coordinate system. 
The grid size is variable and the surface is taken as the 
line between the region of conductivity A, which is a free 
space region, and the region of conductivity B. Anomalous 
conductivity regions which may be of arbitrary shape are 


represented here by the different letters. 


The equation to be solved in all regions for the 


E-polarization case is of the form 


yv2Ex = in2Ex (3.24) 


where n2=4Wvw. 


If Ex = ftig then 


v?ft+ig?g = in?f-n2g (3225) 


and equating real and imaginary parts gives: 


y -f = -n@g (3.26) 


y-g =e ieodae (327) 


If a small region of the mesh is considered as 
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The general two-dimensional model and 


coordinate system. 
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Pievscuaced inerig. 3.3, (3.26) and (2.27)emtstabe 
Satisfied at each point, and in particular at point "0". 
Four conductivities occupy the quadrants surrounding the 
point "0", The mesh sizes vary, and in general 


d1#d24d3#d*. At the point "0" (3.26) and (3.27) become: 


(Vet) 9 o-#(-n2¢)° (3.28) 
and 
(UeO pero (= T),°. (3.29) 


A pair of finite difference equations is 


obtained by making use of Taylor's Theorem which yields 


ESO ety (of) y).0 04 4* (172) (j.<f79 y=) © 
ASE pce Rees aie 0) 


Peet t/a ea Get (2) ie)-f/3 25) ° 


atl? Neer (3.31) 
Come (dy Vio «dst (172) (0) 42 iy < yo 

ids ye ch we, (3.32) 
ee Oe (degre) .0 3 eal 2 alae 7g =) 

(ae es (3. 33) 


and similar equations for g1!,92,93,9%. 


If higher crder terms are neglected, equations 
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Pig. woes Notation used for grid points, dimensions and 
conductivities of the regions surrounding 


poznt fO#% (After Jones and Pascoe, 1971.) 
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(ene9) andy (3529) 4can) be expressedeassa pair of finite 


difference equations: 


£9 {1/ (d?) +1/ (d2) 24+1/(d3) 2+1/(d%) 2} - n2go 


or 


£4. {1/(d*) 24[ 1/ (dt +d3) ]. (1/43. - 1/d41)} + 


pe wee etlelyaldetd* yo) i017 d*s = 17d2) ae 


£3. {1/ (d3) 2+{ 1/(d14+d3) ]. (1/d1 - 1/d3)} + 


£*. {1/ (d*) 2+[ 17 (d2+d*)-}. (1742. -— 174d4)} (3.34) 


re 
£0. (J1/[4 ]2) - n2go = £1p1 + f2p2 + £3D3 + #4D4 (3.35) 


c 


where 


yee = 


ja 


Ds = 


U7 (G4) es Ey (at ees ye}. (17d? 
TAGE ATE Ay (heeds) ois (47/d* 
1/ (43) 24+ 17 (d2443) ]. (1/0! 


Ad etl 17 (detasjy (ide 


1/dt), 
T7aey 
1/d3), 


1/a4). 


A similar equation can be written for the 


go. ()1/ta‘ 32) + n2f0 = gipDi + g2D2 + g3D3 


‘ 


+ 9%D*. 


(3.36) 
(3.37) 
(6728) 


(3.39) 


g values: 


(3.40) 
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Equations (3.35) and (3.40) must be satisfied at each 
interior point of each region. These two equations can be 
solved simultaneously at point "0" for £9 and g® where up- 
to-date values of 5! and gt are obtained from the previous 
iteration. Equations (3.35) and (3.40) must hold for each 
of the surrounding regions. 


That eis * 


£10()1/(.4° 92) ~ (n1)2,g10 = f11p1 + f12p2 

(4 gnapa + guans (3.41) 
e205 ta! J2) - (n2)2,g20 = £21p1 + f2ep2 

(+ £23p3 + fees (3.42) 
£30()1/[d‘ Jz) - (n3)2.g30 = £311 + £32p2 

: + £33p3 + £34p4 (3.43) 
£40(S1/ta° J) (EGCG OU Sas ees seas 

: + £43D3 + £44)4 (3.44) 
g10() 1p a‘ 32) + (n1)2,.£19 = gl1p!1 + gi2p2 

‘4 gtaps + gisps (3.45) 
g20(Jifta! je) + (n2)2.£20 = g2ipt + g22p? 

‘  agesps + g2aps (3.46) 
g30(S1/[a\ J2) + (n3)2.£30 = g3ipt + q32p2 

+ g33D3 + g34*D4 (3.47) 
gso(yatatyey + (nt) 2.240 = gaapi + ge2ne 
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where the first superscript! indicates the conductive 
region considered and the second superscript refers to the 
particular point of interest. The underlined values are 
"fictitious" values and the boundary conditions for the 
interfaces between the regions of different conductivity 
allow these values tc be expressed in terms of known 


values. 


The internal boundary conditions for the E- 
polarization case are that both the tangential and normal 
components of H are continuous across any interface. 
These two components may be expressed in terms of Ex from 


equations (1.7) and (1.8). 


Bye (L/w). (957/92) — (178)25 (9.9/9 Z) (3.49) 
and 
Hz = -(i/w) QfAy) + (1/W) GGA ¥Y)- (3.50) 


The condition for continuity of the tangential components 
applied to each boundary leads to the finite difference 


equations: 
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These equations allow us to express the fictitious values 


of (3.41) to (3.48) in terms of known values. 


Adding 


(3.41), (3.42), (3.43) and (3.44) and making use of the 


fact that 
fers 
grb = g 6 


we obtain 


RCP dle LT tCt + fece VW ES Cee tet 4, 


Similarly, adding (3.45), (3.46), (3.47) 


obtain 


and (3.48) 
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~B£° + Ag® = giC1 + g2C2 + g3c3 + g*c4 (3.55) 
where 

A= 4 ye Ay (dasye (3. 56) 

4 

Cle-pup2 (3.58) 
ce = 4pe (3.59) 
26-42 (3.60) 
C*e=o4D4. (3.61) 


Equations (3.54) and (3.55) are simultaneous equations 


which must be solved for £9 and g®. 


The observed electric field E(obs) can be 
obtained by multiplying Ex by the time dependence exp(iwt) 


and taking the real part of the result. This gives 


Ex(obs) = Re [ (f+ig)exp(iwt) ] = f.cos(wt) 


-g.Sin (wt). (32.02) 


Similarly for the magnetic field components: 
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Hy(obs) = Re [ (i/w).(QEx/)z) ] 
=-(1/w).{@f£/)z).sin(wt) + Qg/)z).cos(wt)}, 


(3.63) 


HZ(obs) = Re [ (-i/w).(gExfy) ] 
=(1/w).{QfA y).sin(wt) + Wg y).cos(wt)}. 
(3.64) 


The phases of these three components may be calculated as 


follows: 


Phase of Ex(obs) = Arctan{{[f.sin(wt) + g.cos(wt) ] 


/(£.cos(wt) - g.sin(wt) ]}, (3.65) 


Phase of Hy(obs) = Arctan{{ (@f/) z).cos (wt) 
- Q9/2z).sSin (wt) J/[- (f/f z) .sin(wt) 


- (99/92) -cos(wt) }}, (3.66) 


Phase of Hz(obs) = Arctan{{- (jf/) y) .cos (wt) 
- QagAy).sin(wt) J/f Qf£/) y)-sin (wt) 


+() g/J) y)-cos(wt) ]}. (3.67) 


3.2.b Combined Solution for Non-Uniform Source and 


Laterally Nen-Uniform Subsurface 
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From the method of Peltier and Hermance (1971), 
the electric field due to a symmetric non-uniform source 
over a horizontally layered earth may be calculated at any 
point (y,Z) where <«<y<+e and z>-j{h]. Values taken from 
such calculations may be used to give the boundary values 
in the Jones and Price (1970) and Jones and Pascoe (1971) 
iterative method of solving for the fields associated with 


lateral inhomogeneities in conductivity. 


When the two methods are combined, the 
conductivity configuration must have identical layering at 
both boundaries for |y]| large, since the Peltier and 
Hermance method applies only for a horizontally layered 
earth. Also, the anomalous conductivity region should be 
approximately three skin depths of the surrounding medium 
from the bottom and side boundaries. This latter point 
has been discussed by Pascoe and Jones (1972). When 
considering the location for the upper boundary of the 
finite difference mesh in the numerical method it is 
desirable to place it at such a height that any 
perturbations due to lateral inhomogeneities are 
negligible and such that it is below any ionospheric 
conductive region. In the work that fcllows the current 
source is chosen at z=-110 km., and the upper boundary of 


the finite difference grid at z=-50 km. 
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The soluticn, after the boundary values have 
been obtained as outlined above, is calculated in the same 


manner aS in the preceding section 3.2.a. 
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4. AN INVESTIGATION OF THE SOURCE EFFECTS IN SYMMETRIC AND 
NON-SYMMETRIC IONOSPHERIC CURRENT DISTRIBUTIONS 


4.1 The Symmetric Current Source Distribution 


4,.1.a Mathematical Formulation of the Boundary Conditions 


The coordinate system which is used throughout 
this section is given in Fig. 4.1. The conductor occupies 
the half-space z>0, so that its boundary is the plane z=0. 
The source current is assumed to flow parallel to the 
surface of the earth at height z=-{h{ and is sinusoidally 
time dependent. The source current is symmetric about the 


line of maximum intensity. 


In section 3.1.b it was shown that for a 
Gaussian source over a horizontally layered earth, the 


elementary electric field in the nth layer can be written 
Ex" = {An.exp(-P.z) + Bn.exp(P.z)}-.G(s).cos (Sy) (4.1) 
where (P)2 = s@tiwu.Vn. 

For a sheet current source with a Gaussian distribution 


function at height z=-[|h{ and a two-layered earth below, 


Peltier and Hermance (1971) show (as given in 3.23) 
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G(s) = -(ilewuk27.s) .exp(sz0-s2k2/2). (G2) 


where Ile is the maximum current intensity, and k is the 


standard deviation cf the source. 


To complete the solution, An and En must be 
determined and the eigenfunction for Ex” given above must 
be summed over all values of s for each point where 


a value for the electric field is desired. That is 


& 
Ex(y;,Z) = f can (2) .exp (-P.2) + Bn(P).exp (P.Z) } 


~G(Ss).cos(sy) ds. (4.3) 


The half-space above the earth is taken as the 
first layer and if A, =1 then G(s) is a frequency 
dependent parameter of the source. By, (P) is thought of as 
a reflection coefficient which represents the contribution 
to the total electric field above the earth from the 
subsurface layers. For a two-layered conductor By =0 so 
that Ex3 will not become infinite with depth. The other 
An and Bn values can be obtained since the tangential 
components of Ex and Hy must be continuous across layer 


boundaries. This leads to 


By amet (Ba1p7 (R41) (4.4) 
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where 

R= Ss. (Q+1)/( P, (Q-1) J (4.5) 
and 

GG= 442, Pa) IP sP5) ) exp {2=P, 4) (4. 6) 


which are identical to equations (3.11), (3.12) and 
(3.13), where d is the thickness of the first conducting 


layer. Also, we will have 


A, = {(1#B,).0}/(0+1), (4.7) 
Be TV AS bow ie (4.8) 
As = {(A, +B, -exp(2.P, d)}.expl (FP, -B).d]. (4.9) 


This general method of course may be extended to 
more layers. However, we will be concerned only with a 


two-layered earth. 


Bete bevariation of (Ex, [Hy |) andy | Hz duestor ay Non— 


Uniform Source above a Layered Conductor 


To understand better the problem of induction by 


a non-uniform source, it is useful to study the variation 
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of the field component profiles with height and frequency 
over a layered earth. A two-layered earth has been 
considered with upper layer of resistivity 100 ohm-m and 
lower layer of resistivity 10 ohm-m. The current source 
with a Gaussian distribution in intensity and standard 
deviation of 240 km. which varies with periods of 10 sec, 


103 sec and 105 sec is placed at height z=-{h]|=-110 km. 


Fig. 4.2A shows the field profiles for the 
amplitudes of the three components Ex, Hy and Hz for 
different heights above the surface of the conductor for 
ten seconds period. This figure clearly illustrates that 
above the surface of the conductor there is increasing 
uniformity of the field components as the surface is 
approached. That is, the source effect diminishes with 
increasing z. This uniformity of the field components 
near the surface is evident at all periods (Figs. 4.2B, 


4.2C). 


As the period increases, the field components 
become more uniform in the y direction, which reflects the 
change in the fundamental wave number of the field. Also, 
as the period increases, the skin depth in the conducting 
regions increases and the effect of the induced currents 


is to make the field more uniform. 
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Pig. 452 {Hz{, |[Hy{ and |Ex{ profiles along the 
surface of the conductor and at heights of 25 
km., 50 km., 75 km. and 100 km. above the 
conductor. (A) 10 "sec... period ,2 (By? 103 
sec. period (C) 105 sec. period. All 
components are normalized with respect toa 
point on the surface of the conductor at 200 
km. from the center of the source. For each 
component the bottom profile is the surface 
profile and the components increase in 
Magnitude with height to the top profile at 


100 km. above the surface. 
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The field component profiles conform well with 
Maxwell's equations. Referring to Fig. 4.2A, the spacing 
between the curves of |Ex|, which are calculated 
for different heights at egual intervals above the 
conductor and so are a measure of J{Ex|/jz, increases with 
decreasing z. Correspondingly, |Hy| increases 
With height. The spacing between the {Ex| curves depends 
on the horizontal distance of the field point from the 
center of the source. The maximum spacing occurs directly 
under the source (at 0.0 km.) and the minimum spacing 
occurs at + 200 km. {Hy{| will then decrease toward + 200 
km. glEx{/AVy steadily increases from a minimum directly 
under the center of the source to maxima at + 200km., and 
so {Hz| increases from a minimum of zero under the center 


of the source to maxima at + 200km. 


It is evident that the field distribution does 
not closely resemble the spatial distribution of the 
current source. The induction effects in the earth are 
principally responsible for this. Price (1965) has placed 
approximate limits on the source field dimensions. He 
estimates that the spatial wavelength of natural fields 
will be from approximately 400 km. to 40,000 km. Since 
the field distribution does not closely follow the spatial 


distribution of the source, this implies that even for 
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rather narrow current sources the spatial wavelength of 
the field still remains within the limits set by Price. 
This 2s evident from Fig. 4.3 where it is shown that for 
current distributions of different widths the electric 
field distribution at z=-100 km. deviates further from a 
Gaussian curve as the source becomes narrower. Also, as 
the period increases the field tends to become more 
uniform. This trend is illustrated in the figure for 
periods of 10 sec. and 1000 sec., and it becomes more 


pronounced for longer periods. 


It is the difference between the field 
distribution and the spatial distribution of the current 
that makes it imperative to compute the distribution of 
the field due to a particular non-uniform current source 
over an appropriately layered subsurface, and to obtain 
from such a calculation the boundary values for insertion 


into a model with lateral discontinuities. 


4U.1.c The Solution for a Laterally Inhcmogeneous Structure 


and Apparent Resistivity Curves. 


After the application of the above method to 
obtain the boundary values for the finite difference 
method, solutions were obtained for a two-dimensional 


model with a lateral inhomogeneity as illustrated in 
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4.3 


E-field profiles at z=-100 km. for current 
sources of standard deviation 240 km., 120 
km. and 60 km. for 10 sec. and 1000 

sec. period compared with corresponding 
Gaussian curves. Normalized with respect to 


the point under the center of the source. 
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Fig. 4.4. The position of the source with respect to the 
inhomogeneous structure is shown in Fig. 4.5. The source 
is located a z=-110 km. and is displaced relative to the 
center of the inhomogeneity by 720 km. The standard 
deviation of the source is 240 km, and so the center of 
the inhomogeneity is 3 standard deviations from the center 
of the source. The field distributions due to the non- 
uniform source over a multiple layered earth without the 
inhomogeneity were calculated at positions corresponding 
to the boundary of the mesh as illustrated in Fig. 4.4. 
These were then inserted into the iterative finite 


difference program and solutions obtained. 


Apparent resistivity curves for different 
positions under the Gaussian source over the layered earth 
without any inhomogeneity are shown in Fig. 4.6 and are 
compared with the apparent resistivity curve for a uniform 
source. It can be seen that at long periods the apparent 
resistivity curves spread away from the uniform source 
curve. This effect was observed by Peltier and Hermance 


(1971) for their models. 


In Fig. 4.7 are given the apparent resistivity 
curves calculated from the Cagniard (1953) formula for the 


conductivity structure shown in Fig. 4.4 for 
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Fig. 


4.4 


Conductivity configuration for the lateral 
discontinuity model. H-values; horizontal 
grid sizes in km.; K-values; vertical grid 
sizes in km. Sigma represents the 
conductivity .in-e. m.u., so that for region 
B, resistivity=0.25 ohm-m; region C, 


resistivity=250 ohm-n. 
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Fig. 4.5 Source location relative to the inhomogeneous 


structure. Distances in km. (Not to scale). 
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Fig. 


4.6 


Apparent resistivity curves for Gaussian 
source over layered earth compared with 
uniform source. The curves are for different 
positions under the source, Q=directly under, 
0.5=one-half a standard deviation from 


the center etc. 
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Fig. 


4.7 


Apparent resistivity curves for different 
positions over the lateral inhcmogeneity. 

(a) left edge of inhomogeneity, Gaussian 
source. 

(b) center of inhomogeneity, Gaussian source. 
(c) right edge of inhomogeneity, Gaussian 
source. 

(d) center of inhomogeneity, uniform source. 


(ce) edges of inhomogeneity, uniform source. 
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three different positions above the inhomogeneity for both 
Gaussian and uniform sources. At shorter periods the 
apparent resistivity varies widely with position for the 
Gaussian source. The apparent resistivity calculation at 
the center of the inhomogeneity tends to approach the 
value of the higher resistivity material. At the edges of 
the inhomogeneity the calculations of apparent resistivity 
give values more closely related to the resistivity of the 
lower resistivity material. For longer periods the non- 
Symmetry of the inducing field above the inhomogeneity 
causes a spreading of the resistivity curves at the edges 
of the inhomogeneity. Apparent resistivity curves for the 
uniform source, as in the layered case, follow the non- 
uniform source curves at short periods but deviate from 


them at longer pericds. 


4.2 The Non-Symmetric Current Source Distribution 


4.2.a Mathematical Formulation of the Boundary Conditions 


Tlilustrated in Fig 4.8 is the coordinate system 
which will be used throughout this derivation. The source 
is a sheet of current which flows parallel to the surface 
of the earth at height z=-{h{|, is sinusoidally time 


dependent and is non-symmetric in spatial distribution 
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with the general intensity distribution of the form 


Y-e€xXp (py) - 


In the two-dimensional E-polarization problem 
with the period sufficiently long so that displacement 
currents may be neglected, Maxwell's equations reduce 


to the diffusion equation as given previously in (1.9). 


By a separation of variables, solutions of the 
form Ex(y,z) = Y(y).Z(z) are obtained. For a horizontally 


layered earth 


Y(y) = C(s,y°)-.cos(sy) + D(s,y®).sin(sy), (4.10) 
and 

Z(z) = An(P).exp(-P.z) + Bn(P).exp(P.2Z) (4.11) 
where (P)2 = s2tiwu.Vn, and s, y® and Vn are the spatial 


wave number, current line source position and the 


homogeneous conductivity of the nth layer respectively. 


The elementary electric field in the nth layer 


can be written as 


EX (y,Z) = {An(P).exp(-P.z) + Bn(P).exp(P.2)} 
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« {C(Ss,y°).cos(sy) + D(s,y°).sin(sy)}. (4.12) 


The elementary magnetic field components 


Hy (y,Z) and Hz (y,zZ) are obtained by taking the curl of 


(4.12). 
Hy (y-2) = -(i/fwu). P{An(P).exp(-P.z) - Bn(P) 

exp (P.z)}. {C(s,y9) .cos(sy) 

teU(Se yo) sin (Sy) } (4.13) 
HZ (y,z) =. +(i/wu). S{An(P).exp(-P.z) + Bn(P) 


-exp (P.z)}.{C(s,y°) .sin(sy) 


— D(S,v¥?). cos (Sy)R.. (4.174) 


If the half space above the earth is defined as 
the first layer, V,=0, Pi =Se and A, (P)=1, then C(s,y) and 
D(s,y°) are parameters which characterize the source and 


are dependent upon the frequency and source position. 


The solution for a current source which has a 
general non-uniform spatial distribution can be obtained 
from the superposition of suitably weighted line 
current source solutions. In order to accomplish this, 
the constants C(s,y9) and D(s,y®) in (4.12) must be 
calculated for a line source. Assuming that the 
conductivity of the whole earth is zero, then B,=0 and the 


horizontal component of the magnetic field Hy! would be 
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given by Ampere's law for a line current source: 
Beer 2/2229) { (2-20) 2s (yy oye) Cer?) 


where z9=-{h{ and y° is the horizontal position of the 


line current source. From (4.13) with V, =0, B,=0 and 


1 
PSs the general sclution for Hy is obtained by summing 
all possible solutions of type (4.15) over the allowable 


Tange of s- (0<s<@): 


w@ 
Hyit = -(iywu). f Ss. {exp (-Sz) }. {C(s, y9) .cos (sy) 
& 
+ D(s,y®).sin(sy)} ds = -(1e/2M) 


elZ=20) 7b (2-22 )et (y= 0) et" (4.16) 


Using the inverse Laplace transforms, 


00 
(ae /{ (2220S ty yc) Ss} = fexp (sz) .exp (295) 
0 0 
~cos(hty=yoes) Eds = [exp (-zs) .exp (20s) 
G 
- {cos (sy°).cos(sy) + sin(sy°).sin(sy)} ds, 


(4.17) 


the integral equation for C(s,y®) and D(s,y®) can be 


solved if we set 


Cis, yo) = C'(s).cos (sy?) (4.18) 
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and 

D(s,y®) = C'(s) .sin (sy °) (4.19) 
where 

C'(s) = -(iwule/2Ws) .exp(z%s). Sea) 


For the non-symmetric current intensity distribution 


Ix(y°) = Ie. (y9-b) .exp {-a(y°—-b) }.u(y9-b) (4.21) 


where Ie, b and a are constants and u(y®-b) is the unit 
step function, we assume that the source is a composite of 
elemental line sources having intensities described by 
(4.21). The magnetic field intensity that such a source 


at (y®,z°) contributes to the observed total field at 


Oe 2 ees 


Hx1(element) = -(1e/2W). (y9-b) .exp {-a (y°-b) } 


ai (Z>aoy/@ (2 2oper (yy oie (ye su). (4.22) 


The total magnetic field Hx1(y,z) can ke obtained by 
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integrating over all y® (b<y%@); 


- 
Bat (YZ) 0= - (102M) « { (yo-b) .exp {-a (yo-b)} 


-{(z-2z°) /[ (2-2) 2+ (y-y°) 2]} dy®. (4.23) 


As before, 


co 
- (i/wu) « fs. exp (-sz) « {V(s) cos (sy) 
& 
+ W(s).Sin(sy)} ds (4.24) 


Hx? (y,Z) 


where V(s) and W(S) are Source constants which can be 
obtained in the same manner as C(S,y9) and D(s,y°). The 
source constants will not be a function of y® (source 
position) since this is the variable over which the 
integration has been taken in (4.24). Equating (4.23) and 


(4.24), we have 


co 


~(i/wu) . f's.exp (-sz). {V(s) -cos (sy) + W(s) 
e - 00 
PSLDASY i. ds: — - (Ze 72m) « f (y°-b) .exp {-a (y°-b) } 
(2-29) (2520 my eS pea yon (4.25) 


If in (4.16) both sides are multiplied by (y®-b).exp{- 


a(y°-b)} and then integrated over all y® we have; 
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-{C(s,y°) .cos(sy) + D(s,y®) .sin(sy)} ds dy? 
(at 
= ER (y°-b) .exp{-a(y-b)}. { (z-z9) 


7 (Ase yiet(y-yoje ly ydy 2: (4.26) 


Comparing (4.25) and (4.26) and solving for V(s) and W(s) 


gives 

V(s) = -(iwule.,exp(sz°) /2Ms) . {{ (a2-s2) «cos (sb) 
~2as.sin(sb) ]/(a2+b2) 2} (ieay) 

and 

W(s) = -(iwule.exp(sz9) /2Ws) . {[ 2as.cos (sb) 


+(a2-s?2).sin(sb) J/ (a2t+s2) 2}. (4.28) 


Summing over ali values of s, the electric field 


in the nth layer is 


N eG 
Ex (y,zZz) = f (an(e) .exp (-P.2) + Bn(P).exp(P.Z) } 
0 


- {V(s).cos(sy) + W(s).sin(sy)} ds. (4.29) 


The magnetic field ccmponents are 


rr) 
Hy “(y,z) = -(4/aa) « f P {an (P) exp (-P.2) - Bn (P) 
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-exp (P.z)}.{V(s) .cos(sy) + W(s).sin(sy)} ds 
(4. 30) 


00 
uz (y,2z) = +(i/wu). f stan(P) .exp(-P. 2) + Bn(P) 
oC 
e€xp(P.z)}.{V(S).Sin(sy) - W(s).cos(sy)} ds. 


(4.31) 


Completion of the solution requires An(P) and 
Bn(P) to be found for each layer. The half space above 
the earth is taken as the first layer, and if A, (P)=1 then 
V(s) and W(s) are frequency dependent parameters of the 
source. These quantities have been previously given in 


section 4.1.a. 


4.2.b Variation of |{Ex|, |Hy]|! and |{Hz{ due to a Non- 


Symmetric Source and Apparent Resistivity Curves. 


It is interesting to compare the field component 
profiles with height and frequency over a layered earth, 
and examples are shown in Fig. 4.9. A two-layered earth 
with the thickness of the upper layer equal to 50 km. and 
resistivity of 100 ohm-m and lower layer of resistivity 10 
ohm-m has been considered. The non-symmetric current 
source has an intensity distribution of (y-b).exp{(- 
1/1b1). (y-b)} with b=-480 km., and is placed at a height 


of z=-|{h{=-110km. The field component profiles for 
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Fig. 


4.9 


{Hzt, {Hy{ and |Ex{ profiles along the 
surface of the conductor and at heights of 25 
km., 50 km., 75 km. and 100 km. above the 
Gonductoir. (A) 10 sec. period, (B) 103 

sec. period (C) 105 sec. period. All 
components are normalized with respect to a 


point on the surface at the origin. 
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periods of 10 sec., 103 sec. ,and 105 sec. areishown in A, 
Band C of Fig. 4.9 respectively. These figures display 
the amplitudes of the three components |Hz{, {Hy{ and {Ex| 
for different heights above the surface of the conductor. 
The maximum amplitude in each case is for the maximum 
height, with amplitude decreasing as the surface is 


approached. 


The non-symmetric current distribution of the 
source is clearly reflected in the field component 
profiles. The field components differ greatly for 
positive and negative y, decreasing rapidly in the region 
of negative y above which the source intensity is zero. 
While the fields are in general non-uniform there 
are similarities to the field components of the symmetric 
Gaussian source given in section 4.1.b. The fields become 
increasingly uniform as the surface is approached. This 
effect is evident at all periods. The source effect 
therefore diminishes with increasing z just as in the 
Gaussian source case. As the period increases the field 
components become more uniform in the y direction, which 
reflects the change in the fundamental wave number of the 


field. This effect was also noted in the Gaussian source 
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Fig. 4.10 for various positions under the source on the 
surface of the conductor. Fig. 4.10A and Fig. 4.10B 

show the apparent resistivity curves in the negative and 
positive y regions respectively at intervals of 200 km. to 
a Maxima of + 1000 km. The curves of Fig. 4.10B resemble 
the Gaussian source curves (see Fig. 2, Peltier and 
Hermance, 1971). However, the curves of Fig 4.10A 
contrast sharply at longer periods. This effect can be 
attributed to the discontinuity of the current source in 


the negative y region. 


4.2.c The Solution fcr a Laterally Inhomogeneous Structure 


and Apparent Resistivity Curves. 


After the application of the above method to 
obtain the boundary values for the finite difference 
method, solutions were obtained for a two-dimensional 
model with a lateral inhomogeneity as illustrated in 
Fig. 4.11... The position of the source with respect to the 
inhomogeneous structure is shown in Fig. 4.12. The source 
is located at z=-110 km. and is displaced relative to the 
center of the inhomogeneity by + 920 km. The intensity 
distribution of the source remains the same as described 
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Fig. 4.10 Apparent resistivity curves for y.exp(py) 
source over layered earth compared with 
uniform source (dashed curve), A-negative y 
region, B-positive y region, V, =100 ohnm-n; 
Vv, =10 ohm-m. 
(a) at origin 
(b) 200 km. from origin 
(c) 400 km. from origin 
(d) 600 km. from origin 
(e) 800 km. frem origin 


(f) 1000 km. frem origin 
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Fig. 
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Conductivity configuration for the lateral 
discontinuity model. H-values; horizontal 
grid sizes in km.; K-values; vertical grid 
Sizes in km. Sigma represents the 
conductivity ine. m. u., so that for region 
E, resistivity=0.25 ohm-m; region C, 


resistivity=250 ohm-m. 
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Fig. 4.12 Source location relative to the inhomogeneous 


structure. Distances in km. (Nct to scale). 
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Apparent resistivity curves for different 
positions under the source over the layered earth without 
any inhomogeneity are shown in Figs. 4.13A and 4.13B. 
Again at long periods the apparent resistivity curves 
spread away from the uniform source curve, and Figs. 4.13A 


and 4.13B differ due to source asymmetry. 


Fig. 4.14 gives the apparent resistivity curves 
for the conductivity structure shown in Fig. 4.11 for 
three different positions above the inhomogeneity for non- 
symmetric and uniform sources. Parts A and B of this 
figure refer to regicns y<0 and y>0 respectively. As in 
the Gaussian source case, for shorter periods the apparent 
resistivity varies widely with position. The apparent 
resistivity value at the center of the inhomogeneity tends 
to approach the value of the higher resistivity material. 
At the edges of the inhomogeneity the calculation of 
apparent resistivity gives values more closely reiated to 
the resistivity of the lower resistivity material. For 
longer periods the non-symmetry of the inducing field 
above the inhomogeneity causes a spreading of the 
resistivity curves at the edges of the inhomogeneity. 
These general trends hold true for both parts A and B. As 
before the positive and negative y region curves differ 


substantially. As in the layered case, the apparent 
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Apparent resistivity curves for the y.exp(py) 
source over layered earth without any 
inhomogeneity compared with uniform source 
(dashed curve). A-negative y region; B- 
positive y region, V, =-25 ohm-m; a oo ohn- 
mM. 

(a) at origin 

(b) 200 km. from origin 

(c) 400 km. from origin 

(d) 600 km. from origin 


(e) 800 km. from origin 


(£) 1000 km. frem origin 
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Fig. 4.14 Apparent resistivity curves for different 
positions over the lateral inhcmogeneity. A- 
negative y region; B-positive y region. 

(a) left edge cf inhomogeneity, y.exp(py) 
source. 

(b) center of inhomogeneity, y.exp (py) 
source. 

(c) right edge of inhomogeneity, y.exp (py) 
source. 

(d) center of inhomogeneity, uniform source. 


(e) edges of inhomogeneity, uniform source. 
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resistivity curves for the uniform source deviate from the 


curves of the non-symmetric source at longer periods. 
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5. ELECTROMAGNETIC INDUCTION IN THE EARTH BY AN APERIODIC 
NON-UNIFORM CURRENT SOURCE 


5.1.a The Two-Dimensional Double Current Systen 


Little attention has been given to the problem 
of electromagnetic induction by a non-uniform source which 
is temporally aperiodic and also exhibits spatial 
variation with time. Price (1950) examined aperiodic 
inducing fields and considered the effects due to a sudden 
increase in the inducing field and the solution for any 
other time variation could then be obtained. The special 
case of a stationary line current parallel to the surface 
of the earth in which a sudden increase in the line 
current intensity tock place was considered by him. Also, 
Wait (1954) considered the transient response of the 
tangential electric field for a unit step increase in the 
magnetic field. The calculations by Wait (1954) which are 
applicable to non-uniform sources, were made for 
subsurface layers of different conductivity and 


susceptibility. 


Kisabeth and Rostoker (1971) show that the 
fields associated with polar magnetic substorms are 
aperiodic and indicate that the structure of polar 


magnetic substorms is controlled by the intensification of 
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the currents at the northern and/or southern borders of 
the auroral electrojet. They propose that quasi-periodic 
intensifications of the northern border of the electrojet 
are responsible for the jagged structure noted on 
magnetograms during the substorms. It is therefore of 
interest to consider the effects at the surface of a 
finitely conducting earth due to multiple current system 


distributions. 


Within the limitations of the two-dimensional 
problem, a double current system may be studied by an 
adaptation of the method of Peltier and Hermance (1971) 
and the results of section 4 for investigating non-uniform 
sources. The two current systems, which may be spatially 
displaced relative tc one another and which may be made to 
vary aperiodically with time, can have maximum intensities 
which may occur at different times and so the current 


system as a whole will vary spatially. 


In the present discussion the two current system 
distributions flow akove an earth with a horizontally 
layered conductivity distribution. One system has a 
Gaussian current intensity distribution, exp(-y2/2k2). 

The other system has a current intensity distribution of 
yeexp(-ay). Each system may vary aperiodically with time 


by summing the solutions for different frequencies through 
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the use of the Fourier transform. The spatially varying 
solution is obtained from the superposition of the 


separate aperiodic time solutions. 


5.1.b Method of Solution 


The coordinate system is shown in Fig. 5.1. The 
Gaussian current source is symmetric about the origin (see 
Fig. 5.2) and the current flows parallel to the surface of 
the earth at a height of 110 km. In the following, this 
source will be referred to as the symmetric source. The 
second source, with current intensity distribution of the 
form y.exp(-ay) will be referred to as the non-symmetric 
source. The current associated with this latter source 
also flows parallel to the surface of the earth at a 
height of 110 km. The non-symmetric source has been 
displaced 700 km. in the positive y direction from the 
symmetric source (see Fig. 5.2). Values for the magnetic 
field components due to the two current systems are 
obtained at intervals of 100km. in the region -1000 


km.<y<1000 km. 


For the E-polarization case, certain symmetric 
and non-symmetric current distributions yield solutions 


for the two-dimensional diffusion equation which results 
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when the period is sufficiently long so that displacement 
currents may be neglected. In the E-polarization case 
(Jones and Price, 1970), this equation is the same as 


CIR Dc 


As shown in section 4 for a horizontally layered 
earth and for an inducing field with sinusoidal time 
dependence, the electric field in the nth layer can be 


written as 


Ex” (Viz) - fone. exp(-P.z) + Bn(P).exp(P.Z) } 


cae cos(sy) + W(s).sin(sy)} ds (5. 1) 
where the constants An(P) and Bn(P) are defined by (4.4) 
through (4.9) and (P)2 = s@tiwu.Vn. 


For a symmetric source over a two-layered earth 


it has been shown in section 4.1.a that 


V(s) = - (ilewuk 2Ws) .exp(sz°-s2k2/2) (5.2) 
and 
W(s) = 0 (oa) 


where Ie is the maximum current intensity, and k is the 
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Standard deviation at the source. Referring to section 


4,.2.a, the non-symmetric current source distribution is 


Bx (Y 20s Bey O by exp {=a (y9=b) 4 ut yo—b) (5.4) 


where b and a are constants, u(y9-b) is the unit step 


function, and 


V(s) = -(iwule.exp(sz°) /2Ws) . {{ (a2-s2} .cos (sb) 
-2as.sin(sb) J/(a2+b2) 2} (5.5) 
and 
W(s) = -(iwule.exp(sz°) /2Ws) . {[ 2as.cos (sb) 
+ (a2-s2) .sin(sb) ]/(a@+s?) 2}. (5.6) 


The magnetic field is obtained by taking the 


curl of (5.1) and the components are expressed as: 
cs a) 
Hy (y,z) = -(i/wu). f P {an (P) .exp (-P.2) =—sb0(P) 
(4) 
eexp (P.Z)}. {V (Ss) .cos(sy) + W(s).Ssin(sy)} ds 


(5.7) 
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N 
Hz’ (y,z) = +(iywu). f s{An(P).exp(-P.z) + Bn(P) 
0 
e€xp (P.z)}.{V(s).sin(sy) - W(s).cos(sy)} ds. 


(5.8) 


The aperiodic solution is obtained for the 


magnetic field values by writing 


co 
Hy (y,t) - fax (y,W) -F(w).exp(iwt) dw, (5.9) 
and 
00 
Hz(y,t) = f Bz (yew) -F UW) -exp (ive) dw (54510) 
GO 


where F(w) is the Fourier transform of the time variation 


€(t) of the current source, 
oO 
F(w) = (172m) feo .exp(iwt) dt. (5.11) 
“0d 


Hy(y,t) and Hz(y,t) are the magnetic field responses as a 


function of time. 


The two magnetic field components of both 


sources were calculated along the surface of a 
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horizontally layered earth with two layers. The upper 
layer has thickness 50 km. and resistivity of 100 ohnm-n. 
The second layer is of infinite depth and resistivity 10 


ohm-m. 


The symmetric source has a current distribution 


I = Ie.exp(-y2/2k2), (5.12) 


where the half-width (k) is 240 km. The non-symmetric 


source has a current distribution 


I = Ie, (y-b) .exp(-a[ y-b ]) .u(y-b) (31:3) 


where b=700 km., a=2.1x10-4 km.-! and u is the unit step 
LUNCELOUM. (Refer to Fig. 5.2 for relative source 


positions.) 


The time variation of both current distributions 


is the same: 


£ (ty) = exp{-d-<. (t—q) =}. (5. 14) 


For 0<t<9 minutes d=0.2222 and for 9<t<63 minutes d=0.040. 


The time variation is shown in Fig 5.3. 
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The Fourier spectrum of the current distribution 
time variation must be band limited to less than 1 Hz or 
the approximation that displacement currents may be 
neglected will not apply. The frequency spectrum of the 


time function chosen is negligible above 0.1 Hz. 


The time variation of the magnetic field 
components Hy and Hz over the region of interest on the 
surface is shown in Fig. 5.4 for the symmetric source and 
Fig. 5.5 for the non-symmetric source. Starting with time 
t=0 a total time span of 63 minutes is shown at time 


intervals of one minute. 


5.1.c The Two-Dimensional Aperiodic Spatially Varying 


Double Current Solution. 


As outlined before, the apericdic spatially 
varying source is obtained by the superposition of the 
symmetric and non-symmetric magnetic field component 
solutions. The non-symmetric source lags the symmetric 
source by ten minutes in time. The time variation for the 
composite source is shown in Fig. 5.6 with a time span of 
74 minutes and time intervals the same as in Figs. 5.4 and 
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The variaticn of Hy (below) and Hz (above) 
Spatial profiles with time for the 
superimposed symmetric and non-symmetric 
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respect to the origin produces an Hy magnetic field 
component at the surface which is also symmetric about the 
Origin. The Hz cross-over point (the point at which Hz 
equals zero) coincides with the origin as can be seen in 
Fig. 5.4. The non-symmetric source, having most of its 
current concentration positioned to the right of the 
region where the calculations of the magnetic field 
components were obtained, shows a high degree of non- 
symmetry in both the Hy and Hz magnetic field components 
as is evidenced in Fig. 5.5. The fact that the current 
concentration is to the right of the region of interest is 
clearly indicated by the position of the maximum value of 
Hy and the cross-over point of Hz. When the two sources 
are superimposed the magnetic field component 
configuration "moves" spatially with time. Hy exhibits 
two maxima in time at different positions along the 
surface as can be seen in Fig. 5.6. The cross-over point 
of Hz shifts to the right with time as can be seen in 
Figs. 5.7, 5.8 and 5.9.1, These are Spatial profiites of the 
Hy and Hz components for intervals of time four minutes 


apace Starting at t=0. 


The non-symmetry in the surface profile of the 
Hy component (Fig. 5.7) for the first two time intervals 


is due to the fact that the contribution of the non- 
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symmetric source was held constant for the first 10 
Minutes. The initial values of Hy at t=0 of Fig 5.5 were 
used for the non-symmetric source contribution of Fig. 5.6 
for t=0 minutes through t=10 minutes. In Fig. 5.6 for 
t>63 minutes, the symmetric source was held constant at 
its final value (t=63 minutes, Fig. 5.4). The same 
procedure was followed for the calculation cf the Hz 
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Kisabeth and Rostoker (1971) and Kisabeth (1972) 
have observed a latitudinal movement of the HZ cross-over 
associated with polar magnetic substorms, and it is 
apparent that this work supports their conclusions with 


reference to the double current sources. 
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6. RECOMMENDATIONS FOR FURTHER RESEARCH 


The three-dimensionality of the physical 
electromagnetic induction problem causes difficulty if the 
direct application of the previously discussed methods of 
sections 4 and 5 to actual observed data is attempted. As 
pointed out by Onwumechilli (1967) allowance must be made 
for return currents, such that the forward current equals 
the return current. This can be accomplished by combining 
several symmetric and non-symmetric sources such that 
magnetic field profiles fit the observed data. 

(e.g. One symmetric source can be allowed to flow parallel 
to the x axis with two non-symmetric sources, one at 
either side, flowing antiparallel.) If two aperiodic 
sources which allow for return currents are then 
superimposed as in section 5, it should be possible to 
make the spatial shift of the magnetic field 


components fit the actual observations. 
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